load data80;
subplot(2,1,1)
plot(cf);
h1=legend('Actual wind speed');
set(h1,'Fontsize',12);
xlim([0,1000]);
subplot(2,1,2)
plot(qx);
h2=legend('One-day-ahead wind speed prediction of NWP');
set(h2,'Fontsize',12);
ylabel('Wind Speed (m/s)');
xlabel('Step (15min)');
xlim([0,1000]);

plot(e);
ylabel('Error (m/s)');
xlabel('Step (15min)');
xlim([0,1000]);
set(gca,'FontName','Times New Roman','FontWeight','bold','Fontsize',12);

subplot(2,1,1)
autocorr(e);
subplot(2,1,2)
parcorr(e);
narxop=[error;output];
save narxop narxop
load narxop
subplot(2,1,1)
plot(narxop(1,:));
h1=legend('Error of NARX model');
set(h1,'Fontsize',12);
xlim([1,1000]);
set(gca,'FontName','Times New Roman','FontWeight','bold','Fontsize',12);
subplot(2,1,2)
plot(narxop(2,:));
h2=legend('Predictive correction');
set(h2,'Fontsize',12);
xlabel('Step (15min)');
ylabel('m/s');
xlim([1,1000]);
set(gca,'FontName','Times New Roman','FontWeight','bold','Fontsize',12);
% load data80test;
% % subplot(2,1,1)
% % plot(cf,'r');
% % h1=legend('Actual wind speed');
% % set(h1,'Fontsize',12);
% % xlim([0,1000]);
% % subplot(2,1,2)
% plot(qx,'r');
% h2=legend('One-day-ahead wind speed prediction of NWP');
% set(h2,'Fontsize',12);
% ylabel('Wind Speed (m/s)');
% xlabel('Step (15min)');
% xlim([0,1000]);